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Abstract 

We study some dynamical properties of a family of two-dimensional cellular 
automata: those that arise from an underlying one dimensional sand automaton 
whose local rule is obtained using a latin square. We identify a simple sand 
automaton T whose local rule is algebraic, and classify this automaton as having 
cquicontinuity points, but not being equicontinuous. We also show it is not 
surjective. We generalise some of these results to a wider class of sand automata. 

1 Introduction 

In |CFM07j , the authors introduce the family of sand automata: these are dynamical 
systems : A^ d — > A^ , where A is an countably infinite alphabet, that satisfy 
certain constraints (see Section [2] for all definitions). In turn, with the appropriate 
topology T put on A zd , (A zd , <£) is topologically conjugate to a cellular automaton 
$q : Sq ->■ Sq where Sq C {0, l} zd+1 is a subshift of finite type. In |DGM09| 
the authors ask if any of the cellular automata arising from sand automata are 
chaotic. This is a particular instance of the more general project of finding chaotic 
higher dimensional cellular automata. In this article we study a family of sand 
automata and their dynamical properties. 

In Section 12.51 we define linear sand automata: these are automata whose local 
rules are built using a group endomorphism <f> : A 2r+1 — > A, where the finite al- 
phabet A is viewed as a cyclic group. We say that local rules are 'built', because, 
unlike cellular automata, the initial configuration has to be relativised before the 
local rule is applied, and the output of the local rule is added to the initial con- 
figuration. In general we work exclusively with one dimensional sand automata (ie 
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those whose corresponding cellular automata act on two dimensional configuration 
space), whose local rules have radius 1. We frequently work with one fixed sand 
automaton, F : A 1, — > A z , whose local rule is built using the rule 7 : Z| — >■ Z5 
defined by 7(2-1 xo x±) = X-\ + x\, which is the local rule for the famous XOR 
cellular automaton (albeit on a different alphabet), and whose dynamical properties 



have been studied extensively, for example in |MM98j . [CFMM97] and [IQN83 



Although we are interested in the dynamical properties of the cellular automaton 
i n practice we often work with <£, since all of the dynamical properties that 
we discuss are topological invariants. In Section [31 we show that T is not surjective 
- recall that a cellular automaton which is not surjective cannot be transitive. In 
Section [37Ll we identify a proper T-invariant subspace, Q, and show that for points x 
in Q, computation of T n {x) is simple, for all n. We identify other sand automata for 
which Q acts similarly - this accounts for about 17% of all radius one sand automata, 
all of which cannot be transitive. 

In Section [J] we show that T is not equicontinuous, by showing the existence of 
(many) vertical inducing points. In fact all radius one sand automata whose local 
rule tables have, for some non-zero m, the value m appear in each column (or each 
row), have vertical inducing points, and so are not equicontinuous. This means 
that at least 99% of all sand automata are not equicontinuous. Despite not being 
equicontinuous, we show that T has equicontinuity points by finding a blocking word 
for r. This completes the classification of T according to the classification scheme 
in |DGM09j . 

In Section I3~2l we generalise the definition of a vertical inducing point to that of a 
local rule constant point. Both vertical inducing and local rule constant points have 
easily computable $ n iterates. We identify a T-invariant subspace G*, in which all 
points are local rule constant. We define a subspace Q' containing Q* , and conjecture 
that Q' is an attractor for T, in that limn^,^ <i(r n x, Q') = whenever x G A z . 



2 Preliminaries 

2.1 Notation 

If A is a countable alphabet, let A + the set of finite concatenations of letters from A; 
elements of A + are called words. If d is a natural number, we let X = A z := {x = 
{(x n ) ngZ d : x n € A for each n € we call X the d-dimensional configuration 

space on the alphabet A, and elements of X are called configurations. If x G A z , 
we sometimes write x = . . . X-i.xqX\ .... If x £ A z , m, n € Z, m < n let X[ mn ] : = 
x m i m+ i . . . x n . We will be working in dimensions 1 and 2, and in dimension 2 only 
when A is finite. 

We let A T ■= {— r, . . . ,r}, and, if A is given, let A = A U {—00,00}. Loosely 
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speaking if d = 2 then x represents a 3-dimensional sandpile where x n represents 
the number of sand grains at "location" n. If x n = oo there is a source at location 
n and if x n = — oo there is a sink at location n. 

Let L C A z<1 be a finite set. If c = {q : c\ E A, for 1 € n L}, define the 
cylinder set [c] = {x : x\ = c\. for 1 € L}. If is finite the Cantor topology - 
generated by the cylinder sets, as L and c vary - makes A^ d compact. If A = Z this 
is no longer the case. Because of this it is convenient to embed Z z in {0, 1}^ 2 , and 
equip it with the subspace topology T that it inherits from {0, l} z2 : this is what is 
done in |DGM09| . Define the map e : Z z -> {0, l} z2 by 



(e(x))ij 



1 */ j < SCi 
i/ j > x { 



In other words, the amount of sand X{ at location % in x is recorded in the ith column 
of e(x), where there are ones for all column entries that are indexed by at most Xj, 
and zeros above entry x%. If Xj = oo, the i-th column of e(x) is constantly 1, and 
similarly if = — oo the i-th column of e(x) is constantly 0. 

The transformation e : X — > e{X) is bijective and its image e(X) is a subshift of 



finite type, where the unique forbidden word is 
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If {0, 1} Z is equipped with the 



Cantor topology then since e(X) is closed in {0, 1} , it is also compact. Let T be 
the topology on J = Z z inherited from the subspace topology on e(X) C {0, l} z2 . 



2.2 Dynamical systems 

We recall some basic definitions. If X is a topological space and F : X — > X be con- 
tinuous, then (X, F) is a topological dynamical system. (X, F) is transitive if for all 
nonempty open sets U, V C X, 3n > such that F~ n (U)nV ^ 0. (X, F) is equicon- 
tinuous at x G X if Ve > 0,35 > 0,Vy € £ 5 (a;),Vn > 0, d(F n (x), F n (y)) < e; and 
(X, F) is equicontinuous if -F is equicontinuous at all points in X. (X, F) is sensitive 
to initial conditions if 3e > 0, Vx € X, V<5 > 0, 3y € ^(x), 3n > 0, d(F n (x) , F n (y)) > 
e. Finally (X, F) is positively expansive if 3e > 0, Vx ^ y,3n > 0,d(F n (x),F n (y)) > 
e. A point x 6 X is periodic if F n (x) = x for some n; if n = 1 then x is /ixed We 
call the smallest such n the period of x. The point x € X is eventually periodic if 
Bp, n £ N such that F p+n (x) = F p (x); we call p the preperiod. A dynamical system 
is ultimately periodic if Vx € X, x is eventually periodic. 

A dynamical system is chaotic if is sensitive, transitive and has a dense set of 
periodic points. In |BBC + 92| it is shown that if (X, F) is transitive and has a dense 
set of periodic points then (X, F) is sensitive. In |CM96] . the authors show that 
if (X, F) is a cellular automaton (see below for definitions) then transitivity alone 
implies sensitivity. 
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2.3 Cellular and sand automata 

Let X = A^ where A is a finite alphabet. A cellular automaton F : X — > X is 
a map defined by a local rule f : A 2r+1 — > A, so that for each n E Z, (F(x)) n = 

n _ r , . . . x n + r ). The simplest example of a cellular automaton is the shift map 
a : X — >■ X, where (<r(x)) n = x n +i. If X = ^4 zd , a cellular automaton can be 
defined similarly: if L r := {1 : 1 = (Zi,^), |4| < is the 2-dimensional box of 
radius r, and / : A Lr — > »4 is a map, then a cellular automaton F : A 7 " 1 —> A 1,2 
is defined as (F(x)) n = /({x n+ i : i E L r }). There are two shift maps on two- 
dimensional lattice spaces X = A^: let (<7ff(x)) n = aj n +(i,o)) (°v( x ))n = £ n +(o,i) 
be the horizontal and vertical shifts respectively. In [Hed69j . Hedlund showed that 
F is a cellular automaton if and only if F is a continuous, shift commuting map. 

Let X = Z z . Like cellular automata, sand automata are defined in [CFM07] as 
functions <E> : X — > X which have a local rule <f>, except that the output of the local 
rule is added to the original entry. First, in [CFM07] the authors define a sequence 
of "measuring instruments" of precision r. If n E N and m G Z then the measuring 
tool of precision r at reference height m is the function /3™ : Z — l A r where 



+oo if n > m + r, 

— oo if n < m — r, (2-1) 
n — m otherwise. 



~2r 

Let ^ : A r — > be given. Define 



<J>(x), 



Xj + 4>{/3^(xi_ r ), . . . ffifa-i), fipixi+i), . . . Ppfa+r)) if Xj G Z 
Xj if Xj = ±oo 



Thus 3>(x)j differs from Xj by at most r + 1, and is a function of the cylinder 
i. Note that 4> is applied not to xu_ rti+r -\ but to a relativised version of this 
word. 

As a book-keeping device we define a projection LT : X — > {A^ r } z by II(x)j = 
^(xi-r), . . . /3^(xj_i), f3'f. l (xi + i), . . . j3p{xi +T ). Clearly LT is not injective. We refer 
to II(x)j as a gradient tuple, and later in this article we use Li, . . . L r to describe the 
gradients to the left of the reference position and R± , . . . R r to describe the gradients 
to the right of the reference position. 

Example 2.1. This is Example 12 in [CFM07], which emulates the behaviour of the 

original model defined in |BTW 88j . Define a 1 dimensional sand automaton F whose 

— 2 

local rule <f> : A\ — » A\ is given by: 



<j>(a,b) 



1 if a = oo, b 7^ — oo 
— 1 if a 7^ oo, b = — oo 
otherwise 
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F has the property that a grain of sand falls to the right (and only the right) if 
the right neighbour is at least 2 smaller. If the number of sand grains in the initial 
configuration x is finite, then F n {x) is eventually fixed. 

Define the vertical map p : X — > X where 



p(x)i 



Xi + 1 if \xi\ < oo, 
Xi if \xA = oo, 



and say that 3> is vertical commuting if = p($(x)). Also, say that <I> is 

infiniteness conserving if 3>(x)j = ±00 -£4> Xi = ±00. Note that all sand automata 
are shift commuting, vertical commuting, and infiniteness conserving; in fact this 
characterises them, as shown in Theorem 17 of [CFM07] : 

Theorem 2.2. $ : X — > X is a sand automaton if and only if $ is continuous, 
shift commuting, vertical commuting and infiniteness conserving. 

The number r in the definitions of cellular and sand automata is called the radius. 
In this article we only consider sand automata of radius 1. 

2.4 Modelling sand automata as cellular automata 

Using the injection e : Z z — > {0, l} z2 , we can transform a sand automaton into a 
2-dimensional cellular automaton as done in [DGM09] , Letting Sn^ denote e(X), 
define $m = e o $ o e _1 : — > Sny With this notation, we have the following 
lemma, whose proof is straightforward: 

Lemma 2.3. $>/i\ commutes with both the vertical and horizontal shifts, and if X 
is endowed with the topology T, then (X,F) = (Sny&ns). 

2.5 Linear sand automata 

If |^4| = 2k + 1, then A can be viewed as the additive group Z, 2 k+i = 
{— k, . . . , 0, . . . , ky. Recall that F : l^k+i ~^ ^ffc+i ^ s a ^near cellular au- 
tomaton if its local rule / : ^k+i ~* ^2fc+i is a group homomorphism. Thus 
f{x.[i-r t i + r]) = a Xi_ r + aiXj_ r+ i + ... + a 2r +iXi+ r where a; and Xi £ Z 2 fc+i for 
< % < 2r + 1. The topological properties of F depend on the coefficients af. for 
example if the a-s are relatively prime, then F is topologically transitive, so that 
many linear cellular automata are chaotic. 

Next we describe how to define a sand automaton using a linear cellular au- 
tomaton. Recall we only consider sand automata of radius 1. In this case we can 
display the local rule in terms of a local rule table. A local rule table has a row for 
each possible left gradient L and a column for each possible right gradient R. Entry 
(L,R) of the table is cp(L,R). For example, Figure [1] is the local rule table for the 
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sand automaton defined in Section [2.31 The local rule table is applied after each 
gradient pair is calculated. For example, if x = . . . , 4, -2, 1, . . . then (H(x))o = (^), 
so L = oo and R = — 1, so that (O(x))o = xq + 1 = 3. 
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Figure 1: The local rule table for the sand automaton 0. 

We now define linear sand automata of radius r. Recall that to define (F{x))q 
first we relativise, sending (x-i, xo, x\) to (f3f°(x-i), f3f°(xi)) 6 «4f. Let / : Ai — > 
Z5 be the bijection defined by 



(/(*))« 



2 
1 




if Xi 
if Xi 
if Xi 
if Xi 
if Xi 



00 

1 





-1 

-00 



Given a group homomorphism cf)* : - 
is a linear sand automaton if the local rule 



5, we say that a sand automaton $ 

(p* o / o n. 

Example 2.4. Let 7* : Z| — >• Z5 be defined by 7*(x, y) = x ®y. 

Let r be the linear sand automaton with local rule 7 = 7, o / o n. Then V is a 
radius 1 linear sand automaton. The local rule table in Figure [2] corresponds to the 
group homomorphism 7*. Note that here the rows and columns are indexed by Z5. 
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Figure 2: The local rule table for V 

In this article we will often be working with V, even though we are interested 
primarily in (see Lemma 12. 3p . For, has radius r = 7 (in fact the local 



6 



neighbourhood can be a 7x3 rectangle) grid, so it is often more practical to work 
with the radius one T. Figure [6] contains the local rule table for Fp-y 

3 Non-surjectivity of T 

A non-surjective sand automaton F has Garden of Eden states: these are configu- 
rations that have no .F-pre-image. A non-surjective cellular automaton cannot be 
chaotic, since it cannot be transitive. The surjectivity of sand automata is shown 
to be undecidable in [CFM07| , In this section we show that the sand automaton 
r is not surjective, and generalise this result to some other one dimensional sand 
automata. 

Recall that a sand automaton <I> : X — > X is surjective on a set Y C X if for 
each y G Y, ^(y') = y for some y' G Y. A configuration is finite if all configuration 
entries are finite, and only finitely many entries are non-zero; let J- denote all such 
points. Similarly, let V denote the set of all cr-periodic configurations, all of whose 
entries are finite. In [CFM07j (Proposition 3.14), the following was shown: 



Lemma 3.1. [fCFMO^ J Let <I> : X — > X be a one- dimensional sand automaton. 
Then 

1. is surjective on V if and only if ^ is surjective, and 

2. If <!> is surjective on F, then $ is surjective. 

We show that T is not surjective by first showing that there is a word which has 
no predecessor word under T. In the following proof the notation [n]s denotes the 
projection of n G Z to Z5. 

Lemma 3.2. Let w = (100, 3, 2, 100). Then there is no word y such that T(y) = w. 

Proof. Suppose that y = yo Hi 2/2 V3 2/4 2/5 is a word of length 6 such that T(y) = w. 
Then 

98 < yi < 102, 1 < y 2 < 5, < y 3 < 4, and 98 < y 4 < 102 . 

This implies that 0*( yi ) = (3f(y 4 ) = 00, ie /(/?f ( yi )) = /(/3f (y 4 )) = /(°o) = 2. 
Next under the action of V we have: 

V2 + f(J3?(yi)) + fW! 2 (2/3)) = 2/2 + 2 + /(/3f (2/3)) = 3 (3.2) 

2/3 + /(/3f (2/4)) + fW! s (2/2)) = 2/3 + 2 + /(/?f (2/2)) = 2 (3.3) 
Projecting (|3.2[) and (|3.3[) into Z5, we have 

Ms e 2 [y 3 - y 2 } 5 = -2 (3.4) 
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Ns © 2 [y 2 - y 3 } 5 = 2 (3.5) 
Adding (HOD and (03]) we get 

Ms © Ms = i (3.6) 

Thus the only possibilities for [2/2)5 and [7/3)5 are given by: 
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This implies that the only possibilities for 2/2)2/3 are: 
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Each of these cases implies a contradiction to Equation f|3.6|) . 

□ 

Proposition 3.3. T is not V-surjective, so that T is not surjective. 

Proof. Let x be the periodic configuration x = 100, 3, 2, 100, TOO, 3, 2, 100. Lemma 
13.21 implies that there does not exist y such that T(y) = x. Lemma I3T1 implies that 
T is not surjective. □ 

3.1 Surjective subsets 

While r is not surjective, in this section we identify a proper, closed, T-invariant 
subspace, Q. We then identify sand automata $ for which Q is also ^-invariant. Let 

Q := {x : \xi — Xj_i| > 2, Vi} . 

We shall show that each member of Q has a T-predecessor (though not necessarily 
in Q). We will explain in geometric terms how V acts on G, identify subsets of Q 



8 



that are a possible attractor for T, and generalize these properties to other sand 
automata. 

It is clear that / o Ii{Q) c {( 2 ), (" 2 2 ), (_ 2 2 ), {zl)} Z , and that / o U(G) is the 
subshift of finite type whose transition graph G is shown in Figure Re-label 
(2) = -l, (- 2 ) = 0", (_ 2 2 ) = + , (Z 2 2 ) = 1, and let Y G denote the image of / o U(Q) 
under this labelling. Let {—1,0, 1} Z be the full shift on three letters and let p : 
{(T,0 + , 1,-1} -> {-1,0,1} by defined by p(0 - ) = p(0+) = 0, p(l) = 1, p(-l) = 
— 1; then p is a radius local rule for the cellular automaton P : Yq — > { — 1, 0, 1} Z . 
With this notation, the following lemma is straightforward. 




Figure 3: This graph G defines a subshift on its bi-infinite paths. 



Lemma 3.4. // g 6 Q then Ely € Yq such that T(g) = g + P(y)- 



L\R 


-2 


2 


-2 
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-1 



Table 1: The local rule table for T restricted to Q 

On Q the local rule table for V can be compressed to Table [TJ Call 3-tuples such 
that /oII(a, b, c) = (~ 2 ) peaks (centred at b) and similarly /o]I(a, b, c) = ( 2 ) valleys 
(centred at b). Also label 3-tuples such that / o U(a, b, c) = ("g 2 ) up-slopes (centred 
at b) and / o n(a, b, c) = (\) down-slopes (centred at b). With this labelling of 
gradient tuples as geographical features, the action of r n is easily described, once we 
have knowledge of how T acts. 

Proposition 3.5. If g € Q and y € Yq are such that T(g) = g + y , then r™(g) = 
g + ny for all n in N. 

Proof. Here we claim that all of the geographical features are preserved under T. 
If we show this then the proposition follows. First note that according to Table HJ 
|(r(#))n — 9n\ < 1 for each n, whenever g G Q. Let (5-1, go, gi) = (a,b,c); we show 
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that geographical features are preserved at F(g)o; the cases at F(g) n for n ^ are 
identical. 

1. If (a, b, c) is a valley centred at b, then (T(<?))o = go — 1 (see Tabled]). Since for 
n = ±1, (T(g)) n — g n is at least -1, then the valley at b is mapped to another 
valley centred at (T(g))o. The case where (a,b,c) is a peak is similar. 

2. If (a, 6, c) is a down-slope centred at b, then we either have a peak or a down- 
slope centred at a. Thus (T(g))_i > Similarly there is either a valley 
or a down-slope centred at c, so (T(g))i < g\. This means that a down-slope 
centred at go is mapped to a down-slope centred at (r(g))o- The case where 
(a, 6, c) is an up-slope is similar. 

□ 

Note that the proof of Lemma 13.51 shows that j(G) C Q- each geographical 
feature under the action of the sand automaton can only become more pronounced 
or stay the same. However V : Q — > Q is not surjective. Consider g = 0,3, -0,3, 

so that r(g) = g H 1,1, -1,1. If r _1 (g) contained an element g' in Q, then 

g = g' H — 1,1,-— 1,1 by Lemma 13.51 So g' = 1,2, -1,2 Q, a contradiction. The 
next lemma tells us that although T is not surjective on Q, all configurations in 
P(Yg) are used when determining T(Q). 

Lemma 3.6. If y G P{Xg) then 3g G <5 suc/i i/iai T(g) = g + y. 

Proo/. First find an element y* G {©, (" 2 2 ), („ 2 2 ), (l^)} 2 with P(y*) = y. Then 
choose go arbitrarily and follow the instructions given by y* to specify <?_i and g\ 
- for example, if j/q = (-o) t nen choose g\ < go — 2 and > yo + 2. Continue 
this process, using the gradients at locations -1 and 1 to specify g_2 and gi, moving 
outwards from the central cell. The result follows by induction. □ 

The proofs in this section for the sand automaton T relied on the preservation 
of certain geographical features when restricted to the set of configurations Q. We 
now show that knowledge of the entries a, /3, 5 and A noted in Tabled is sufficient 
to extend these results. 

Theorem 3.7. Let the radius one sand automaton $ have local rule table as in Table 
Ul Then $ is peak preserving if and only if a > (3,5, A and $ is valley preserving if 
and only if X < f3,5,a. 

Proof. Assume that $ is the sand automaton described above. We use the notation 
x [n-i,n+i] = ( a ,b,c) and ($(x))r ri _i ) „ + i| = (a*,b*,c*). We prove the statement 
concerning peak preservation - the proof of the statement for valleys is similar. 
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Table 2: The positions in the local rule table that are used to create peak/valley 
preserving sand automata. 

Suppose that a > (3,5, A. Let (a, b, c) be a peak centred at b. Then b* = b + a. 
(see Table [2]). Since a is the largest increase when restricted to Q, then a*, c* < 
a + a, c + a. Thus peaks are mapped to peaks. Conversely suppose that $ is peak 
preserving, and suppose that a < f3. Then the configuration g = . . . — 2, • 2, . . ., 
that has a peak at the central position and an up-slope at the —1-st position, does not 
map the central peak to another peak. This is due to the fact that (<I>(x))o = 2 + a 
and (<J>(x))_i = + j3, where a < /3. Therefore (<J?(x)) - (<&(x))_i > -2. Similar 
configurations can be constructed when a < 5 or a < A using a valley or a down- 
slope at the 1 position. 

□ 

For the results of Lemma 13.51 to hold, which would also imply that Q is $ in- 
variant, it must be both peak and valley preserving as well as being up-slope and 
down-slope preserving. This implies that we are assuming a > (3, 5 > A. The next 
lemma demonstrates that peak and valley preservation implies up-slope/ down-slope 
preservation. 

Proposition 3.8. Let & be a radius sand automaton with local rule table as in Table 
[H If & is both peak and valley preserving then it is also up-slope and down-slope 
preserving, so that Q is $ invariant and not surjective. 

Proof. Let $ be both valley and peak preserving. Then a > /3, 5 > A. We show that 
down-slopes are mapped to down-slopes, the up-slope case being similar. Suppose 
there is a down-slope centred at x n = b, where X[ n _ l n+1 ] = (a, b, c). Then ($(x)) n = 
b + {3. We either have a peak or a down-slope centred at a. Thus(<I>(x)) n _i > a + f3. 
Similarly there is either a valley or a down-slope centred at c. Therefore ($(x)) n _|_i < 
c + (3. Thus a down-slope centred at b is mapped to a down-slope centred at the 
image of b. □ 

If we count the number of sand automata that are both peak and valley pre- 
serving by fixing both a and A and then counting all 5, (3 that do not violate the 
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conditions a > /3, 5 > A, then in total there are 105 radius one peak and valley 
preserving automata. There are in total 5 4 = 625 possible choices for a, /3, 6, A; thus 
the set of sand automata which is invariant on the set Q represents = 0.168 of 
the total radius 1 sand automata. 

4 Equicontinuity and points of equicontinuity 

In this section we investigate the equicontinuity of radius one sand automata. 
In |CFM07j and [DGM09], the authors classify one dimensional sand automata 
as: either sensitive, or nonsensitive without an equicontinuity point, or non- 
equicontinuous with an equicontinuity point, or finally equicontinuous. 

4.1 Vertical inducing points and equicontinuity 

Given that a sand automaton is topologically conjugate to a 2-dimensional cellular 
automaton, we use the following result: 

Theorem 4.1 (Proposition 3.14, [DGM09] ). <I> is equicontinuous if and only if & 
is ultimately periodic, if and only if Vx £ X, <&(x) is eventually periodic. 

In order to classify T we introduce the following definitions. Let n £ N. A 
configuration x is a vertical inducing point of order n for a sand automaton $ if 
<&(x) = p n (x). For example, a fixed point for <E> is also a vertical inducing point of 
order 0. 

Lemma 4.2. If x is a vertical inducing point of order n, then $ m (x) = p mn (x), for 
each m € N. Also e(x) £ SVi\ satisfies ^Jl-,(x) = <7y n (x), for each m € N. 

Proof. If m = 1, then $ x (x) = p n (x) by definition. Suppose that for m = k, 
$ fc (x) = p kn {x). Then 

$ fc+1 (x) = $($ fc (x)) = $(p fcn (x)) = p fc "($(x)) = p kn {p n {x)) = p( fc+1 )"(x). 

Also, for each m, $"K (e(x)) = e($ m (x)) = e(p mn (x)) = a™ n (x.) . □ 
\o) 

Corollary 4.3. If a sand automaton $ admits a vertical inducing point of nonzero 
order, then <I>^ is not equicontinuous. 

Proof. Let n € N. By Lemma 14.21 if x is a vertical inducing point of order n ^ 

then <J>7K (e(x)) = (Ty Tl (x), Vm € N. The only way that e(x) is eventually periodic 
Ko) 

is if x consists entirely of sinks or sources; however in this case x would be vertical 
inducing of order 0. Theorem 14.21 now implies the result; the fact that equicontinuity 
is a topological property implies the second statement. □ 
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In the next theorem we identify a class of sand automata that have vertical 
inducing points. 

Theorem 4.4. Let <3? be a radius one sand automaton. If there exists a nonzero 
m such that m appears in each row, or each column, of <I> 's local rule table, then <1> 
admits a vertical inducing point of nonzero order. 

Proof. Suppose that the nonzero value m occurs in every row of the local rule table 
for the sand automaton Note that if / o H(x) = ((jj*))ieZ> then Lj+i = —Ri 
for each i. Conversely, if ((^ 4 ))j e z G /(n(-X~)), is such that each Li+x = —Ri, then 
there is some x € X such that f(U(x)) = ((^)); e z G f(U(X)). To find a vertical 
inducing point then, it is sufficient to find a cycle (^), • • • , (^) where L i+ i = —Ri, 
for 1 < i < k, L\ = —Rk and such that </>(^ i ) = m, for each 1 < i < k. First 
select any (^) such that 4>{Li,R\) = m. If Ri = —L\ then we are done. If not, 
select (^ 2 ) such that Li = —R\ and cp(L2, R2) = m - this this can be done because 
the local rule table has the value m in every row. If R2 = —L\ or R2 = —L2 then 
we are done; if not continue this process until we find (^) . . . (^) such that for 
some 1 < j < k, Lj = —Rk - this has to occur since the gradient pair entries come 
from a finite alphabet. The desired cycle is . . . (^)- The proof when there is a 
nonzero m in every column is similar except that newly selected gradient pairs will 
come before the current gradient pairs. □ 

Recall that an n x n latin square is a table filled using an alphabet A of size n, 
where each row and each column has exactly one occurrence of each member of A. 
Theorem 14.41 and Corollary 14.31 then imply the following: 

Corollary 4.5. // a radius one sand automaton $ has a Latin square local rule table 
then $ is not equicontinuous. 

Let $ be a radius 1 linear sand automaton. Then the local rule for <1>, denoted 
by (j), is a homomorphism on a cyclic group with group operation addition ©. In 
particular this implies that (j){L, R) = aL © f3R, where a and f3 6 Z5. This implies 
that <l>'s local rule table is either a latin square, or has a nonzero element on the 
anti-diagonal. 

Proposition 4.6. Let § be a radius 1 linear sand automaton. Then <3? admits a 
vertical inducing point, and so is not equicontinuous. 

Proof. It is sufficient to show that there is a word . . . in (Zjj)+ such that 
Li + \ = —Ri for i = 1, . . . k — 1, R n = —Lq and there exists a nonzero m such 
that for each i, j(Li,Ri) = m. Suppose first that there is a nonzero value on the 
anti-diagonal of the local rule table for <3?. Then there is a gradient pair (_ a a ) that 
returns a nonzero value m under 6. In this case, the word ( a ) is the desired one. 
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If all anti-diagonal rule entries are 0, then for each L £ Z5, aL © (3{—L) = 0. This 
implies that aL = (3L, so that a = f3. Thus <j> = aL(BaR = a(L®R). But the local 
rule L(& R corresponds to the linear sand automaton T which has a Latin square as 
a local rule table. Therefore a(L©i?)'s rule table is a permutation of this local rule 
table and so also a Latin square. Corollary 14.51 now yields the result. □ 

Note that linearity is not used in the case where there is a nonzero value on the 
anti-diagonal of the local rule table. Thus we can make the following statement for 
general radius one sand automata: 

Proposition 4.7. // a radius one sand automaton <I> has a nonzero value on the 
anti-diagonal of its local rule table then <3? has a vertical inducing point. 

There are 5 25 radius one sand automata, and since there are 5 20 local rule tables 
with only zero entries on the anti-diagonal, at least 99.968% of the total number of 
radius 1 sand automata have a vertical inducing point, and so are not equicontinuous. 

An example of a vertical inducing point for T is x = 0, 1, 3, 1, • 0, 1, 3, 1, 0. 
The_n_/(n(x)) = (?),(-/), (3 .(i), (J) -©.(I 1 ), (3 ,(i),©, so that r(x) = 
x+1-1 7^ x. In fact the sand automaton T has an infinite number of vertical inducing 
points for each —2 < n < 2. We discuss this in the following section. First though 
we show that while T is not equicontinuous, it does have equicontinuity points, 
putting it in Category (3) of the classification in [DGM09] . We say that a word w 
is blocking for a sand automaton <I> if there 3k, s £ N such that < k + s < |w| and 
such that whenever x and y £ [w]i, then ($ n (x))[ i+M+ | w |_ s ] = ($ n (y))[i+fc,i+|w|- s ] 
for all natural n. 

Proposition 4.8. Let w = 03 230. Then w is blocking for the sand automaton T. 
Thus r has equicontinuity points. 

Proof. First note that the first statement will imply the second. Note also that 
/ o n(w) = (\°), (I 2 ), (I), (Zl), (I)- This block satisfies 7 ((l 2 ), (}), (:£)) = 
(2, 2, 2). Suppose that x, y £ [w]_2- Then (r(x))[_ 22 ] = (x'_ 2 , 5, 4, 5, x' 2 ) and 
(r(y))r_2 2] = W-21 5; 4, 5, y' 2 ), where x'^y^ < 2, for i £ {—2,2}. Under the action 
of r the central three cells increase by 2 and the right most and left most cells can 
increase by at most 2. This implies that (/ o II(r(x)))r_ lj i] and (/ o n(r(y)))[„ 11 ] 
equal (/ o n(xr_ 1) u)), and V thus adds 2 to these three positions. An inductive 
argument completes the proof. □ 

There are 24 nontrivial homomorphisms / : Z 2 — > Z5. If f(L, R) = aL © f3R we 
use the notation (a, (3) to represent /. Eight of these maps have blocking words of 
the type described in Proposition 14.81 These are (a, a), where a / 0, and (—1,2), 
(1,-2), (2,-1), (-2,1). 
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4.2 Local-rule-constant configurations 

The concept of a vertical inducing point is a special case of a more general type of 
configuration x where ($ n (x)) n is easily computable. Let us say that sand automa- 
ton is local-rule- constant at x (or x is a $ local-rule-constant point) if for some 
y G X, $ n (x) = x + ny for each n G N. For example, if x is a vertical inducing 
point, then it is a $ local-rule-constant configuration, with y constant. However the 
set of local-rule-constant configurations is much larger; we now describe a family of 
these points for V. 

Let w = W\ . . . Wfr = (^) . . . (^) G (Z5) 4 " satisfy the properties that Ri = 
—Li + i for i = 1, ... k — 1, and 7*(t0j) is constant i = 1, . . . k; then we say that w is a 
cyc/e segment. In Tabled we list some cycle segments for T. Consider the directed 
graph T~L whose vertices are the cycle segments for V listed in Table [3l and such that 
there is an edge from vertex V to vertex V if and only if the following are satisfied: 

1. If the cycle segment corresponding to V ends with a gradient pair (*) then the 
cycle segment corresponding to V starts with a gradient pair (~ a ). 

2. If the cycle segment corresponding to V ends with a gradient pair (*) and has 
order j, then the order of the cycle segment corresponding V must have order 
at least j. 

3. If the cycle segment corresponding to V ends with a gradient pair and 
has order j, then the cycle segment corresponding to V must have order at 
most j. 

Note that vertices K, L and M in "H are isolated. Define Q* to be the set of all 
configurations x in X such that /(II(x)) corresponds to an infinite path in T~L. In 
this context the infinite loops at K, L and M correspond to T-fixed points in X. 

Theorem 4.9. 7/x 6 Q* , then x is T local-rule-constant. 

Proof. Choose an x G X such that / o II (x) is represented by an infinite path 
V = . . . V-2 V-i ■ Vo V\ V2 ... G H. Let y be the point in X obtained by applying 
7* to the representative in f(U(X)) of V. We claim that T n (x) = x + ny. Suppose 
that Vj corresponds to the cycle segment (/(LI(x)))j j _ 1+ i, . . . (/(II(x)))j j , and it ends 
with a gradient pair of the form (* 2 *). Then in x, a^.+i — Xi j > 2. Geographically 
speaking there is a "steep hill" to the right of X{.. By condition (2), Vi can only 
be followed by an Vi+i whose corresponding cycle segment has order at least that 
of Vi's. Therefore in T(x) the "steep hill", if it changes, can only get steeper. Thus 
/(n(r(x)))j j = /(n(x))j j . Similarly if Vi ends with a gradient pair of the form ( *g) , 
Condition (iii) guarantees that if /(II(x))j. = —2 then /(II(r(x)))j. = —2. This fact 
is true for all j. Finally if V corresponds to the infinite loop at K, L or M, then 
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/(II(x)) is constant and y = • 0, so that T(x) = x in which case T is (trivially) 
local-rule-constant. Thus foIl(Tx.) is also represented by V. By induction it follows 
that x is r-local-rule-constant. □ 



For example, suppose that we want a configuration x that under the action 
of r we have r n (x) = x + ny where y = 2, -1, 1, 1, 1, 2. Then to build such a 
configuration using the cycles from the vertical inducing points we can let /oll(x) = 

©.©.cs.ffi.ca- (-J . (o) > (?) , c 1 ) . cm cm a- ™ s ^ to 

an infinite number of configurations in the pre-image set (/ o IT(x)) . One such 
point is x = 2, 2, 4, 3, 4, -2, 1, 1, 2, 5, 4, 5, 3, 3. 



Label 


Cycle segment 


Order 


A 


CrKXX^) 


1 


B 


(3 


1 


C 


(2)6) 


2 


D 


(30 G) 


2 


E 


C 2 2 ) 





F 


(" 2 2 ) 





G 


r 2 )(-° 2 ) 


-2 


H 


©CD© 


-2 


I 


C 2 ) Co 1 ) (_°i) (i) 


-1 


J 


© 


-1 


K 


(3 





L 


C\) 





M 









Table 3: All of the cycle segments that can be used to create points r n (x) = x + ny. 

Note that the set Q defined in Section 13.11 is contained in Q* . Note also that 
Q* is closed and T-invariant. An interesting question is whether Q* is an at- 
tractor set for the sand automaton T. We have conducted simulations where 
the space time diagrams of several initial configurations are generated. Empir- 
ically what seems to be happening is that the iterates of the initial configura- 
tion converge "almost everywhere" to a configuration in Q* . We describe what 
we mean by this: define the words O = C2) (0) (2) (-2) ' P = Ci) (0) (1) (-2) and 
Q = ©Co 2 ) (-2) ©'^ = (DCo^Ci)©- These sets of words are Aperiodic in 
the sense that if / o n($ n (x) [M+3] ) = O then / o n($ n+1 (x) [M+3] ) = P and if 
/ o n($ ri (x)[j j j + 3]) = Q then / o n($ n+1 (x)[j j + 3]) = R. If we include these in a new 
graph Ti' , then this seems to describes the asymptotic behaviour of more accu- 
rately. The graph %' is presented in Figure[5l This leads to the following conjecture. 
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Figure 4: The graph for the set of all points made up of constant blocks 
which will add a consistent y at each time step. See Table [3] for definitions of 
A,B,C,D,E,F,G,H,I,J,K,L,M. 

Similar to our definition of Q*, let 

Q' : = {x : /(II(x)) corresponds to an infinite path in H'}. 

Conjecture 4.10. The set Q' is an attractor for V, in that if x 6 X, then 
\ixn n ^ 00 d{T n (x),g')=Q. 
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Figure 5: A potential attractor set for T. 
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Figure 6: The local rule table for r7i\. The number on the left of each rectangle 

represents the number of ones in the reference column. An 'x' represents a cell that 

20 

can have any value. The ■ represents the central cell. Only configurations where the 
central cell changes are listed. 



